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ABSTRACT 

We present a model to self-consistently describe the joint evolution of starburst galax- 
ies and the galactic wind resulting from this evolution. This model will eventually be 
used to provide a subgrid treatment of galactic outflows in cosmological simulations of 
galaxy formation and the evolution of the intergalactic medium (IGM). We combine 
the population synthesis code Starburst99 with a semi- analytical model of galactic 
outflows and a model for the distribution and abundances of chemical elements inside 
the outflows. Starting with a galaxy mass, formation redshift, and adopting a particu- 
lar form for the star formation rate, we describe the evolution of the stellar populations 
in the galaxy, the evolution of the metallicity and chemical composition of the inter- 
stellar medium (ISM), the propagation of the galactic wind, and the metal-enrichment 
of the intergalactic medium. The model takes into account the full energetics of the 
supernovae and stellar winds and their impact on the propagation of the galactic wind, 
the depletion of the ISM by the galactic wind and its impact on the subsequent evo- 
lution of the galaxy, as well as the evolving distributions and abundances of metals 
in the galactic wind. In this paper, we study the properties of the model, by varying 
the mass of the galaxy, the star formation rate, and the efficiency of star formation. 
Our main results are the following: (1) For a given star formation efficiency a more 
extended period of active star formation tends to produce a galactic wind that reaches 
a larger extent. If /* is sufficiently large, the energy deposited by the stars completely 
expels the ISM. Eventually, the ISM is being replenished by mass loss from supernovae 
and stellar winds. (2) For galaxies with masses above I^^^^Mq, the material ejected 
in the IGM always falls back onto the galaxy. Hence lower-mass galaxies are the ones 
responsible for enriching the IGM. (3) Stellar winds play a minor role in the dynam- 
ical evolution of the galactic wind, because their energy input is small compared to 
supernovae. However, they contribute significantly to the chemical composition of the 
galactic wind. We conclude that the history of the ISM enrichment plays a determi- 
nant role in the chemical composition and extent of the galactic wind, and therefore 
its ability to enrich the IGM. 

Key words: cosmology: theory — galaxies: evolution — intergalactic medium — 
ISM: abundances — stars: winds, outflows — supernovae: general. 



1 INTRODUCTION 

Galactic winds and outflows are the primary mechanism 
by which galaxies deposit energy and metal-enriched gas 
into the intergalactic medium (IGM)13 This can greatly af- 
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^ Some authors make a distinction between galactic winds, which 
are generated over most of the lifetime of the galaxy and inject 
energy and metals at a steady rate, and galactic outflows, which 
result from violent processes like starbursts, are short-lived, and 



feet the evolution of the IGM, and the subsequent forma- 
tion of other generations of galaxies. Feedback by galactic 
outflows can provide an explanation for the observed high 
mass-to-light ratio of dwarf galaxies and the abundance of 
dwarf galaxies in the Local Group, and can solve various 
problems with galaxy formation models, such as the over- 
eject material at large enough distances into the IGM to eventu- 
ally reach other galaxies. In this paper, we use one or the other 
to designate any material that is ejected from the galaxy and 
deposited into the IGM. 
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cooling and angular momentum problems fsee lBensonll201(il 
and references therein). Galactic outflows can explain the 
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and provide observational tests that can constrain theoreti- 
cal models of galaxy evolution. Local examples of spectacu- 
lar outflows in dwarf starburst galaxies include those of the 
extre mely metal-poo r I Zw 18 (Fequignot 2008i:|jamet et al.l 
I2OIOI ) and NGC1569 jWestmoauette et al ll2009l) . More mas- 



sive spirals, such as NGC7213 jHameed et al]|200ll ). also 
show ev idence of global outfl ows. For a review of the sub- 
ject, see IVeilleux et all (|2005l '). 



1.1 Galactic Outflow Models 

Large-scale cosmological simulations have become a major 
tool in the study of galaxy formation and the evolution of the 
IGM at cosmological scales. These simulations start at high 
redshift with a primordial mixture of dark and baryonic mat- 
ter, and a spectrum of primordial density perturbations. The 
algorithm simulates the evolution of the system by solving 
the equations of gravity, hydrodynamics, and (sometimes) 
radiative transfer. Adding the effect of galactic outflows in 
these simulations poses a major practical problem. In one 
hand, the computational volume must be sufHciently large 
to contain a "fair" sample of the universe, typically sev- 
eral tens of Megaparsecs. On the other hand, the physical 
processes responsible for generating the outflows take place 
inside galaxies, at scales of kiloparsecs or less. This repre- 
sents at the very minimum 4 orders of magnitude in length 
and 12 orders of magnitude in mass, which is beyond the 
capability of current computers. Since we cannot simulate 
both large and small scales simultaneously, the usual solu- 
tion consists of simulating the larger scales and using a sub- 
grid physics treatment for the smaller scale. Cosmological 
simulations can predict the location of the galaxies that will 
produce the outflow, but cannot resolve the inner structure 
of these galaxies with sufBcient resolution to simulate the 
actual generation of the outflow. Instead, the algorithm will 
use a prescription to describe the propagation of the outflow 
and its effect on the surrounding material. 

One possible approach consists of depositing momen- 
tum or thermal energy "by hand" into the system, to 
simulate the effe ct of galactic outflows on the surround- 
ing m aterial (IScannapieco et al.l [2OO1I: iTheuns et al] 



20021: 



20051 : 



Springel & Hernquisi^ '2003'; ' 
Qppenhe imcr & Dave 2006; KoUmcior et al. 200eJ). The 
algorithm determines the location of the galaxies producing 
the outflows and calculates the amount of momentum 
or thermal energy deposited into the IGM based on the 
galaxy properties (mass, formation redshift . . .). Then, in 
particle-based algorithms like smoothed particle hydrody- 
namics (SFH), this momentum or energy is deposited on 
the nearby particles, while in grid-based algorithms it is 
deposited on the neighboring grid points. This will result 
in the formation and expansion of a cavity around each 
galaxy, which is properly simulated by the algorithm. 

A second approach consists of combining the numer- 



ical simulation with an analytical model for the outflows. 
Tegmark et al. (1993) have developed an analytical model 
to describe the propagation of galactic outflows in an ex- 
panding universe. In this model, a certain amount of en- 
ergy is released into the interstellar medium (ISM) by 
supernovae (SNe) during an initial starburst. This en- 
ergy drives the expansion of a spherical shell that prop- 
agates into the surrounding IGM, until it reaches pres- 
sure equilibrium. This model, or variations of it, has been 
used extensiv ely to study the effect of galactic outflows 



on the IGM (iFurlanetto fc L ocb 2001^ : iMadau et al l 
Scannapieco &: BroadhurstI 2001.: Scannapieco et aU 



2001 



scannapieco taroadnursti /uui.: scannapieco et al.l 120021 : 
Scannapieco fc Ohll2004l : iLevine fc Gnedinll2005l; iFieri et al 



2007 , hereafter FMG07; ISamui et aU I2OO8I : iGermain et al 
200£ ; iFinsonneault et al.ll201ol ). In this approach, the evo- 



lution of the IGM and the propagation of the outflow are 
calculated separately, but not independently as they can in- 
fluence one another. The presence of density inhomogeneities 
in the IGM can affect the propagation of the outflow, while 
energy and metals carried by the outflow can modify the 
evolution of the IGM. 

There are several limitations with this second approach. 
In particular, it assumes that the initial starburst, which 
occurred during the formation of the galaxy, is the only 
source of energy driving the expansion of the outflow. First, 
the starburst lasts for a short period of time, typically 
50 million years. Hence, we would expect to observe very 
few galaxies having an outflow. Furthermore these galaxies 
would be just forming and therefore would have complex and 
chaotic structures. Observations show instead that outflows 
are ubi quitous and often or iginate from well-relaxed galax- 
ies fsee lVeilleux et al]|2005l and references therein). Second, 
even though the injection rate of energy is maximum dur- 
ing the initial starburst, the total amount of energy which 
is injected afterward by all generations of SNe and stellar 
winds could be comparable or even more important. Even 
if this energy is injected slowly over a large period of time, 
the cumulative effect could be significant. Indeed, an initial 
outflow caused by a starburst could be followed by a steady 
galactic wind that would last up to the present. The role of 
stellar winds has been mostly ignored in analytical models 
and numerical simulations of galactic outflows. Third, there 
is the possibility that accretion of intergalactic gas onto the 
galaxy might trigger a second starbur st. Finally, a recent 
study JSharp fc Bland-Hawtho"rnll2010l ) suggests that there 
is a signiflcant time delay between the initial starburst and 
the onset of the outflow, something not considered by cur- 
rent models. 

Another important issue is the amount of metals con- 
tained in the outflow, the spatial distribution of metals in 
the outflow, and the relative abundances of the various ele- 
ments. The metallicity of the outflow depends on the metal- 
licity of the ISM at the time of the starburst. The met- 
als contained in the ISM at that time can have several ori- 
gins: (1) metals already present in the gas when the galaxy 
formed, (2) the SNe produced during the starburst, and (3) 
the stellar winds generated by massive stars and AGB ob- 
jects. Hence, the composition of the outflows will depend on 
the epoch of formation of the galaxy (which determines the 
initial metal abundances), as well as the relative amount of 
metals injected into the ISM by Type la SNe, core-collapse 
SNe (Types lb, Ic and II), and winds, and the timing of these 
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various processes. As for the amount of metals injected in the 
IGM, models often assume that it is proportional to the mass 
of the galaxy, and do not provide a description of the distri- 
bution of metals in the outflow and the relative a bundances 
of the elements l|Scannapieco fc Broadhurstll200ll : PMG07). 



1.2 Objectives 

Our goal is to develop a new galaxy evolution model to im- 
prove the treatment of galactic winds in cosmological sim- 
ulations. This model will describe not only an initial star- 
burst and its resulting outflow, but the entire subsequent 
evolution of a galaxy up to the present. It will take into 
account the progressive injection of energy by SNe and stel- 
lar winds (which could cause a steady galactic wind that 
would follow the outflow and last up to the present), and the 
time-evolution of the metallicity and composition of the ISM 
that would directly affect the composition of the galactic 
wind. It will also provide a description of the metal content, 
metal distribution, and chemical composition of the galactic 
wind. This emphasis on the structure and composition of 
the galactic wind is what distinguish our model from recent 
semi-analytical models of galaxy formation, which tend to 
focus on reproducing the properties of the galaxies them- 
selves (luminosity and mass functions, halo properties, disk 
sizes, . . .). For a revi ew of the various semi-analytical mod- 
els, see lBaughl (|2006l ). Our approach combines a population 
synthesis algorithm to describe the stellar content of the 
galaxy, an analytical model for the expansion of the galac- 
tic wind, and a new model for the distribution of elements 
inside the galactic wind. 

The paper is organized as follows: in §2, we describe 
the method we use to calculate the mass and chemical com- 
position of stellar winds and SN ejecta, a key ingredient of 
our algorithm. In §3, we describe the basic equations for the 
evolution of the ISM. In §4, we describe our galaxy wind 
model. Results are presented in §5, and conclusions in §6. 



2 POPULATION SYNTHESIS MODEL 
2.1 Starburst99 



Starburst99 |Leitherer et al.l 19921: Leitherer fc HeckmanI 



1 19951 : iLeitherer et al.1 1 19991 : I Vazquez fc Leithereij l2005l l is 
the population synthesis code we selected to simulate the 
properties of the stellar winds and SNe produced by stellar 
populations with various characteristics. In this project, all 
simulations use an instantaneous SFR to produce a stellar 
population of mass i\fpop ~ 10^ Mq , with a standard IMF as 
defined by [Kro upa ( 2001il . This population is evolved up to 
a final time tflnai = IGyr, with a timestep At = 10*yr (small 
enough to reproduce each evolutionary phase). The evolu- 
tionary tracks are a combination of the Geneva tracks at 
young ages (^ lO^Myr) with the Padova tracks at old ages, 
to better reproduce the mass loss phases of the massive stars 
as well as those of smaller-mass objects. We consider four 
values of the initial metallicity {Zi = 0.001, 0.004, 0.008, and 
solar metallicity 0.02). We also specify that all stars with an 
initial mass Mi of 8 to 120 Af© will end as SNe. Binary sys- 
tems which may also produce SNe are not included in this 
code and not considered here. 




log t [yr] 

Figure 1. Energy and mass ejected by a single stellar population, 
versus time. The SNe mechanical energy (LsNci top panel), stellar 
wind mechanical energy (LsWi middle panel), and mass loss rate 
by stellar winds (Mgw, bottom panel) have been obtained for a 
10® Mq stellar population, using Starburst99. Different metallic- 
ities have been considered as indicated in the top panel. 



The mass we chose for the stellar population, Afpop — 
10^ Mq, is large enough to provide a good sampling of the 
IMF. Yet, it is much smaller than the stellar mass of even 
dwarf galaxies. This enables us to simulate any galaxy SFR 
we want. To simulate an instantaneous starburst with a stel- 
lar mass larger than Afpop, we simply rescale the results pro- 
duced by Starburst99. To simulate a SFR that is extended 
over a finite period of time, we offset the simulations, such 
that during each timestep At, the correct stellar mass is 
formed. In ^S.ll below. we describe the various SFRs consid- 
ered in this study. 

Starburst99 provides, as a function of time, the mechan- 
ical energy produced by SNe and stellar winds, the mass loss 
by stellar wind, and the chemical composition of the stellar 
winds, which takes into account the following elements: H, 
He, C, N, O, Mg, Si, S, and Fe. The code can also predict 
the mass loss by SNe, but for this project we prefer to use a 
more detailed treatment which also gives the chemical abun- 
dances of the SN yields for a more extensive list of elements. 
This treatment is described in the next section. Therefore, 
the ejected mass calculated by Starburst99 will only include 
the contribution of stellar winds. In Figure [T] we plot the 
luminosity (rate of mechanical energy injection) associated 
to SNe, LsNo, and stellar winds, Lsw, and the mass loss 
by stellar winds, Msw, versus time, for stellar populations 
with different metallicities. The initial metallicity has very 
little effect on the SNe luminosity, but greatly impacts the 
luminosity and mass loss by stellar winds. This implies that 
stellar populations that form later will have more powerful 
winds, since they form out of an ISM that has already been 
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enriched in metals by earlier populations. Figure [T] shows a 
sudden increase in wind power around 1 - 2 x lO^yr, just 
prior to the arrival of the first SNe, which is caused by the 
evolved stages of OB stars, i.e. the Wolf-Rayet stars. The 
increase in mass loss after lO^yr is caused by the low-mass 
stars on the AGB. 



2.2 SN Abundances 

Except in situations where the metallicity of the stars is 
larger or equal to solar, th e enrichment of the I SM is domi- 
nated by SNe (see Fig. 8 of lLeitherer et al.lll992f) . Thus, it is 
important to know the composition of the material ejected 
by SNe as a function of the initial mass and metallicity of 
the stellar progenitors. Several groups have used models of 
stellar evolution and nucleosynthesis to produce tables of 
chemical abundances for SN ejecta. These studies are sum- 
marized in Table [1] The first column gives the authors of 
the papers; the second and third columns list the metallic- 
ities and initial masses considered, respectively; the fourth 
column indicates whether or not mass loss prior to the SN 
explosion was included in the model. 

Although very useful, these studies are not fully 
satisfactory for o ur wor k. The first three s tudies 
dWooslev fc WeaveJ 1 19951 : IChieffi fc Limongil l2004l : 
iNomoto et al.ll2006l . hereafter N06) consider a wide range of 
metallicities, from no metallicity (Z = 0) to solar metallicity 
{Zq = 0.02), but are limi ted to initial masses Mj ^ 40 M©. 
The ne xt three studies JWooslev fc Hegeij 2007, he r eafter 
WH07; iLimongi fc Chieffil l2007l : iHeger fc WooslevI |2010| . 
hereafter HWIO) consider initial masses up to 100 - 120 Mq, 
but only one value of the metallicity. Also, three of these 
studies do not include mass loss prior to the explosion, 
which is critical in our models since we include stellar wind 
effects. The study that comes the closest to our needs is the 
one of NOG. This is the only study that covers metallicities 
in the range Z = - Zq and includes mass loss. However, its 
mass range only goes from 13 to 40 Mq . These tables need 
to be extrapolated both at the low-mass and high-mass 
end to cover the full range of SNe progenitor masses. In 
this paper, we assume an IMF with lower and upper mass 
limits of 1 and 120 Mq , respectively, and SNe progenitor 
masses in the range Ali — 8 - 120 Mq . The lower limit of 
8Mq is the most commonly used, but that value is actually 
quite uncertain, and could have an important effect on the 
results (see Fig. [5] below). 

The top panel of Figure [2] shows the mass ejected by 
SNe versus initial mass and metallicity, calculated with the 
tables of N06 (solid curves). An interesting feature is the 
linear relation between Moj and Mi at Z = 0. Also, all the 
curves converge together at Mi = 13 Mq , suggesting that 
for initial masses Mi < 13 Mq, the ejected mass is indepen- 
dent of metallicity. To extrapolate down to Mi — 8M©, we 
used the composition for Mi = 13 Mq, modulated using the 
linear reg ression at Z = calculated with the software Slope 
l|lsobe et a l. 1990). The result is also shown in Figure [2] 

Extrapolation up to 120 Mq can be tricky, especially if 
mass loss is considered. Several approaches have been sug- 
gested in the literature, espe cially at the high-mass end. 
iMartinez-Serrano e t al. (2008) extrapolate the mass ejected 
Mej and chem ical composition of the eject a linearly up to 
Mi = 100 MM. IOppenheimer fc Dwi (|2008l ') assume similar 




Figure 2. Top panel: Total mass ejected by one SN versus the 
initial mass of the progenitor, according to the tables of N06 
(solid lines), WH07 (lower dotted line), and HWIO (upper dotted 
line). Colors represent various metallicities, as indicated. Bottom 
panel: Total mass ejected by one SN versus the initial mass of 
the progenitor after extrapolating down to Mi = 8 Mq and up 
to Mi = 120 Mq, for metallicities Z = (solid black line) and 
Z = Zq (solid green line), after extrapolating over all metallici- 
ties available in N06 (blue and red lines), and after interpolating 
to the metallicity Z = 0.008 used by Starburst99 (purple line). 



ejecta for stars in the range Mi — 3 5 - 100 Mq and also for 
star s in the range Mi = 10 - 1 3 Mm . iTornatore et all (|2007l l 
and IScannapieco et al] (l2005t ) assume that stars with ini- 
tial masses larger than 40 Mq end up directly in a black 
hole, with no ejecta. Using linear interpolation up to Mi = 
120 Mq would be risky. The relation between Mcj and Mi 
becomes metallicity-dependent for Mi > 25 - 30 Mq , im- 
plying that only two or three points would be available to 
extrapolate over a mass interval three times wider than the 
one of NOG. Furthermore at Z = Zq, M^j eventually de- 
creases with increasing initial mass and a linear extrapo- 
lation would eventually lead negative values before reach- 
ing Mi — 120 Mq. Setting the ejected mass to zero for 
Ali > 40 Mq is not a good solution either. Even if the SN re- 
sults in the formation of a black hole, that black hole is never 
as massive as the progenitor fWooslev et al."2 00^; WH07; 
iLimongi & ChicfH 2008; Zlian ct al. 2008; Bclc zvnski et al.l 
|2010|) . Using a copy of the M = 40 Mq at larger masses 
is not realistic either since the tendency clearly shows that 
there is more material ejected for larger initial masses, ex- 
cept at Zq. 

Since none of the three extrapolation methods consid- 
ered so far is satisfactory for extrapolating up to 120 Mq, 
we developed an alternative method. We started by extrap- 
olating the table of NOG at Z = and also Zq, by com- 
bining them with the tables of HWIO and WH07, respec- 
tively. The top panel of Figure [2] shows the Zq relation of 
WH07 and the Z = relation of HWIO (dotted curves). The 
Z = Zq relations of NOG and WH07 differ significantly in 
the range Mi = 20 - 40 Mq . This is likely the consequence 
of using a different prescription for the mass loss during 
the pre-SN phase. In the model of WH07, the ejected mass 
reaches a maximum around Ali = 21Mq, then decreases 
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Table 1. SN yields available in the literature. 



Paper 


Z 


M, [Mq] 


M 


Wooslev & Weaver (1995) 


0, 2 X 10-*^, 2 X 10"*, 0.002, 0.02 


11, 12, 13, 15, 18, 19, 20, 25, 30, 35, 40 


no 


Chieffi & Limonei (2004) 


0, 10"'', 10-"', 0.001, 0.006, 0.02 


13, 15, 20, 25, 30, 35 


no 


Nomoto et al. (2006 ~) 


0, 0.001, 0.004, 0.02 


13, 15, 18, 20, 25, 30, 35, 40 


yes 


Wooslev & Heser (2007) 


0.02 


32 Mi's between 12 and 120 


yes 


limongi & Chieffi (2007) 


0.02 


15 Mi's between 10 and 120 


yes 


Heeer & Wooslev (2010) 





120 Mi's between 10 and 100 


no 



with increasing initial mass, and finally reaches a plateau at 
Mej = 4 Mq , which represents the minimum mass a SN can 
eject. Assuming that this limit is correct, we extrapolated 
the Zq model of N06 until it reaches the model of WH07 
(around Mi — 70 Mq). We then switched to the model of 
WH07 to complete the extrapolation up to Mi = 120 Mq. 
For Z = 0, the results of N06 and HWIO are essentially 
identical in the range A'h = 13 - SOA^q, and very similar 
in the range Mi = 30 - 40Mq. To combine them, we used 
the tables of N06 in the range Mi = 8 - 40 Mq and the ta- 
bles of HWIO in the range Mi = 55 - 120 Mq, Between 40 
and 55 Mq , we interpolated between these two masses. The 
results are shown in the bottom panel of Figure (2] 

At the two intermediate metallicities Z = 0.001 and 
0.004, the values of the ejected masses should lie between the 
values for Z = and Zq (Fig.[2ll. If the mass loss due to stel- 
lar winds were neglected, the masses ejected by SNe should 
be the same as for models with Z = 0. The mass loss associ- 
ated with massive stars is proportional to Z°' , where a ~ 0.6 
- 0.8 ( Vink ot al. 2001; Vink & do Kotor 2005; Krticka 2006; 
iMokiem et al.ll2007l ). We then interpolated to get the mass 
ejected by SNe at any metallicity: 

Mej(M„Z) = Mej(M„0) -r(M„Z)M, M = AZ"^ , (1) 

where A is a constant, and T{Mi, Z) is the lifetime of a star 
of initial mass Mi and metallicity Z. We made the approxi- 
mation that the lifetime does not vary much with metallicity, 
r{Mi,Z) Si r(Mi). Equation ([1} becomes: 

Moj(M„Z) = Mej(M„0)-B(MO^". (2) 

where the function B{Mi) was calculated from the masses 
ejected at Z = Zq: 

BiM.) = M.,iM.,ZQ)-M.,iM.,0) 
^0 

We choose the value a — 0.625, which minimizes the dis- 
continuity at M = 40 Mq . Figure O shows all the models 
of N06 after extrapolation. For Z = 0.004 and 0.001, the 
mass loss is not as extreme as in the case Z = Zq. For that 
reason, we used the composition of the 40 Mq models mod- 
ulated according to the value of Mej for all the initial masses 
Mi > 40 Mq. 

Next we determined the composition of the ejected 
masses, for Z — and Zq. Table [2] shows the hydrogen 
mass, helium mass, and total mass ejected for the WH07 
models. The cases presented are extreme. The mass loss by 
stellar winds is sufficiently large to expel the entire hydrogen 
envelop and most of the helium as well. This indicates that 
these stars have evolved through a Wolf-Rayet stage, result- 
ing in SNe of Type lb or Ic. Hence, most of the material 
ejected is composed of metals. For this reason, we treated 



Table 2. Hydrogen mass, helium mass, and total mass ejected 
by SNe versus initial mass at solar metallicity, according to the 
table of WH07. 



Mi[MQ] 


Mh[Mq] 


MhoIA/q] 


aUIMq] 


40 


0.223 


1.640 


9.74 


50 


5.51 X 10"** 


0.068 


7.94 


60 


1.77 X 10"^ 


0.120 


5.65 


70 


1.16 X 10-^ 


0.152 


4.35 


80 


1.04 X 10"^ 


0.156 


4.34 


100 


1.77 X 10"^ 


0.181 


3.96 


120 


1.39 X 10"^ 


0.172 


3.92 



hydrogen and helium as we treated the total mass ejected: we 
followed the tables of N06 until it reaches the model of WH07 
and then switched to the model of WH07 (Fig. The re- 
minder of the ejected mass was assumed to have the chemical 
composition corresponding to the Mi = 40 Mq model. For 
Z = 0, the complete composition of the Mi — 40 Mq model 
was used. 

Starburst99 requires tables of SNe at metallicities Z = 
0.001, 0.004, 0.008, and 0.02. At this point we have tables 
at Z = 0.004 and 0.02. For Z = 0.008, we interpolated 
between the tables at Z = 0.004 and 0.02 using log Z as an 
interpolation variable. We now have a full set of SN tables 
that provides the ejected mass and composition of the ejecta, 
for all progenitors in the range Mi = 8 - 120 Mq and Z = - 
Zq. Figure [2] shows the ejected mass versus progenitor mass. 

At each timestep during the simulation, Starburst99 
provides the number of SNe and t heir lum inosity (adopt- 
ing that each SN produces lO^^ergs: lMcKeQ,1990. '). The SN 
tables are then used to calculate at every timestep the mass 
and composition of the material deposited into the ISM. For 
simulations with initial metallicity .Z = 0, we use the lumi- 
nosity at Z = 0.001 which is the smallest value considered 
by Starburst99. Figure |4] shows the mass loss rate due to 
SNe for a 10® Mq instantaneous burst. This completes the 
results shown in Figure [T] In the top panel of Figure [5] we 
plot the total mass ejected by SNe versus metallicity, for 
the same population. The dotted curve shows our extrapo- 
lation method, while the solid curve s show the three other 
methods. The "constant" method of lOppenheimer fc Pavel 
(200^ (using the ejecta for Mi — 40Mq at higher masses) is 
the one which resembles our method the most, but there are 
significant differences. This justifies a posteriori the method 
we have developed in this section. 

We assume a minimum progenitor mass of A/i = 8Mq , 
but we also generated other tables using different values. 
The bottom panel of Figure [5] shows the ejected mass for 
minimum progenitor masses of Mi = 6Mq , Mi — 8Mq (the 
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Figure 3. Hydrogen (top) and helium mass (bottom) ejected by 
one SN versus the initial mass of the progenitor, according to the 
tables of N06 (solid lines), after extrapolating down to Mi = 8Mq 
and up to Mi = 120Mq , for Z = Zq . The dotted line illustrates 
the Z = Zq model of WH07. Colors have the same meaning as 
in Figure 12] 
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Figure 4. Mass loss rate by SNe versus time. Number of stars and 
SN explosions are from a Starburst99 model for an instantaneous 
burst with a standard IMF and a total mass of 10** Mq. Mass 
ejected by each SN comes from our extended tables. Different 
metallicities have been considered as indicated. 



value assumed in this paper), and Mi — IOMq. The difTer- 
ences are of order 50%, with the ejected mass being larger 
for a smaller minimum progenitor mass. This shows that 
the particular choice of minimum progenitor mass can have 
a impact on the results, and we intend to investigate this 
issue in more details in future work. Here we focus on the 
case of a minimum prog enitor mass of Mi = 8M q. 

In a recent paper, iHoriuchi et al.l (|201ll ) point to a 
serious "supernova rate problem": the measured cosmic 
massive-star formation rate predicts a rate of core-collapse 
supernovae about twice as large as the observed rate, at least 
for redshifts between and 1, where surveys are thought to 
be quite complete. Several explanations are proposed to ex- 
plain this major discrepancy, including a large fraction of 
unusually faint (intrinsically or dust-attenuated), and thus 
unaccounted for, core-collapse SNe and a possible overesti- 
mate of the star formation rate based on the current estima- 
tors. If indeed this supernova rate problem is real, the SFR 
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Figure 5. Mass ejected by SNe versus metallicity. Top panel: 
results obtained with our extrapolation method for the SN 
mass ejected as a function of the progenitor initial mass, us- 
ing the same stellar population model as in Figure |4] (dotted 
curve), compared with results obtained with different extrapo- 
lation method s (solid curves): from top to bottom: linear ex- 
trapolation ofjMartine^^erran^e^anijSo^ constant approx- 
imation of [O^^enheimerfc el ll2008l'l. and no-ejecta approx i- 
mation of lTornatore et ah I ||2007^ and IScannapieco et al.l ll2005tl . 
Bottom panel results obtained with our extrapolation method, for 
various values of the minimum initial mass. From top to bottom. 
Mi = 6M0, 8Mq, and IOMq. 



(see §5.1) might have to be scaled accordingly. However, all 
simulations presented in this paper start at redshift z = 15, 
and terminate at redshifts between 6 and 9. It is not clear 
that there is a supernova rate problem at these redshifts. 



3 EVOLUTION OF THE INTERSTELLAR 
MEDIUM 

During the evolution of a galaxy, the ISM is constantly en- 
riched by ejecta from stellar winds and SNe. Hence, every 
generation of stars provides an environment richer in met- 
als for the future generations. The level of this enrichment 
depends on the SFR, since the metal production increases 
with the number of stars formed. To simulate this process, 
we designed an algorithm that combines the outputs of Star- 
burst99 with the SNe tables of N06. 



3.1 Initial Conditions 

We consider a galaxy with a total mass Mgai. We assume 
that the ratio of baryons to dark matter in the galaxy is 
equal to the universal ratio, which is a valid assumption for 
the initial stages of the galaxy. The baryonic mass of the 
galaxy is then given by 

Mt = ^A4,, , (4) 

"0 

where Qo and Q,bo are the total and baryon density pa- 
rameters, respectively. We assume that each galaxy starts 
up with a primordial composition of hydrogen and helium 
{X = 0.755 and Y = 0.245). The formation of the galaxy 
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results in a starburst, during which a fraction /* of the bary- 
onic mass is converted into stars. The total mass in stars at 
the end of the starburst is therefore 



M. = /.Mi,. 



(5) 



We refer to the parameter /« as the star formation efficiency. 

We consider three different types of star formation rate: 
instantaneous, constant, and exponential. With an instan- 
taneous SFR, all stars form at i = 0. In the other cases, the 
stellar mass formed M* and the star formation rate M« are 
related by 



M,{t)dt = M^ 



(6) 



where t/ is the final time of the simulation. For a constant 
SFR, we have 



f ' 

"•burst 

where fburst is the duration of the starburst. We usually 
choose a value for M, and solve equation for tburst- For 
an exponential SFR, we have 



M. = 



17 



-t/tc 



(8) 



where tc — 5 x lO^yrs is the characteristic time. Since star 
formation never ends in this case, the stellar mass formed 
depends on the final time tf. Equation ((Sj is valid in the 
limit tf ^ tc. 

There are some caveats about equations Q and ([5]). 
The infalling gas must cool and form molecular clouds which 
then fragment into stars. If star formation is delayed until all 
the gas has cooled, then equation Q would be valid, but it 
is more likely that star formation will start while a fraction 
of the gas has not cooled yet. SNe and stellar winds from 
that first generation of stars will inhibit the formation of 
subsequent generations of stars, by reheating the ISM and 
possibly expelling some of it in the form of a galactic out- 
fiow. Observations of galaxies in the redshift range < 2 < 4 
show that star fo rmation is a very ine fficient process, with 
Af,/Afgai < 0.03 (|Behroozi et al.ll2010l ). Hence, the value of 
Afi, in equation Q is an upper limit, which does not take into 
account the gas lost by galactic outflows. As for the reheat- 
ing of the gas and suppression of inflow, it is implicitly taken 
into account in equation ((5]) by introducing a star formation 
efficiency /*. The mass of gas that was reheated and pre- 
vented from forming stars is Afrohcat = (l~/*)Afi,. In this pa- 
per, we consider star formation efficiencies /« = 0.1, 0.2, 0.5, 
and 1.0. Equations @ and (O then give A/*/A/gai = 0.016, 
0.032, 0.081 and 0.162, respectively. The first two values are 
consistent with observations. The last two values are quite 
extreme, and we considered them mostly to investigate the 
behavior of the model under extreme conditions. 

In this paper, we impose various forms for the SFR to 
investigate the effect of the SFR on the properties of the 
galactic outflows. The effect of feedback and reheating of 
the ISM is all contained implicitly in the value of /*. To 
provide a proper treatment of the aforementioned feedback 
processes, we intend to modify the model such that the SFR 
will be recalculated at every time step, from the physical 
conditions of the ISM gas at that time, taking into account 
the effect of all previous generation of stars. This will pro- 
vide a consistent treatment of feedback and self-regulating 



star formation, and eliminate the need to specify a priori 
a star formation efficiency /, . This will be the subject of a 
forthcoming paper. 



3.2 Evolution of the Mass and Composition of the 
ISM 

Our algorithm tracks the evolution of A/ismxi the mass of 
element X contained in the ISM. At the beginning of each 
timestep, the total mass of the ISM is given by 



A^ISM(t) =^AflSMx(i): 



(9) 



where the sum is over all elements included in the algo- 
rithm, that is all elements from hydrogen (X = 1) to gallium 
(X = 31). These quantities are initialized at the beginning 
of the simulation, and updated during each timestep At, as 
follows. First, we calculate the mass ejected by stellar winds 
and SNe. We include the contribution from all stars present 
at that time, taking into account their current ages, initial 
metallicities, and initial masses: 

MswAt) = J2^'^sw^^n,Zk,Mk)At, (10) 
MsNcAt) = ^MlS!^{n,Z^,Mk)At. (11) 

where Tk, Zk, and Mk are the current age, initial metallicity, 
and mass of population k, respectively. The sums are over 
all the stellar populations that have already formed by time 
t. The superscript SB99 indicates quantities calculated by 
Starburst99. We also calculate the total luminosity produced 
by SNe and stellar winds: 



(12) 
(13) 



Lswit) = Lsw^jrk, Zk, Mk) , 

k 

LsNcit) = -^Inc^ {Tk, Zk, Mk ) . 

k 

We then remove from the ISM the total mass of the 
stars born during that timestep, and the mass removed by 
the galactic wind, and add the material ejected by SNe and 
stellar winds: 

A^isMx {t + At) = AfisMx it) - ^m,sm{!) 

+MGw(t)] At + Afswx (t) + MsNex (i) , (14) 

where Mcwit) is the rate of mass loss by galactic wind (the 
calculation of A/gw is presented in the next section). The 
mass remove from the ISM by the star formation process 
is used to generate several new stellar populations. Each 
population is given an initial mass Mk , an initial metallicity 
Zk equal to the metallicity Z{t) of the ISM at that time, 
and we set the age Tk of these new populations to zero. We 
then recompute the ISM metallicity: 



Z{t + At) = 



A'/isM {t + At) - A/ismh [t + At) 



-AfisMHe(i + At) 



MisM{t + At). (15) 



This expression accounts for the material removed from the 
ISM by the star formation process and the galactic wind, 
and it also considers the enriched material added by stars 
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already formed. Finally, we update the age of every stellar 
population: 



for all k . 



(16) 



These operations are repeated at every timestep in the sim- 
ulation. 

In our model, Mismx is a function of time only. This 
assumes that metals ejected into the ISM by SNe and stel- 
lar winds are instantaneously mixed. This commonly-used 
approximation has the advantage of being simple to imple- 
ment. However, in reality, it will take some finite time before 
the metals are fully mixed. For this reason, the efficiency of 
metal-enrichment of the IGM in our model should be con- 
sidered as an upper limit. 

Starburst99 cannot calculate directly the mass loss and 
luminosities appearing in equations pU|l to p3|) for any ini- 
tial metallicity Zi. It is limited to the values Zi = 0.001, 
0.004, 0.008, and 0.02. We therefore need to interpolate the 
results of Starburst99 in order to get the quantities M^^®^ 
and L^^^^. For values of Zi in the range [0.001,0.02], we 
calculate M^^^® and L^^^^ at the metallicities immediately 
before and after Zi, and interpolate between them, using 
equations of the form 



logM^ 
logL^ 



AlogZ^ + B, 
C\ogZ, + D. 



(17) 
(18) 



In the range [0,0.001], we treat stellar winds and SNe dif- 
ferently. We assume that, in that range, the mass loss by 
stellar wind is proportional to Z" , with a ~ 0.625, as in 
equations ([1]) to ((3|. Hence, 



MlT[Zi^ 
rSBgg/ 



= MS 



\Zi) 



rSB99 



^(0.001) ( 

(oil) 



-V 



0.001 y 

2 X 0.625 



(19) 
(20) 



For SNe, we use the values L^^{Zi = 0.001) at lower 
metallicities. This is a valid approximation because the de- 
pendence of the stellar lifetimes on metallicity is weak, as 
Figure [1] shows. 



3.3 Mass Loss by Galactic Wind 

The presence of a galactic wind enables a fraction of the 
ISM to escape the galaxy and enrich the surrounding IGM. 
Since the wind is generated by the thermal energy deposited 
in the ISM by stars, we expect the mechanical energy of 
the galactic wind to be proportional to the rate of energy 
injection by SNe and stellar winds: 



iMGw(i)V<iw (X L(t) , 



(21) 



where Mqw and Vgw are the mass loss rate by galactic wind 
and the velocity of the wind, respectively. The galactic wind 
will create a cavity expanding into the IGM (see § 4 below) . 
The expansion of the cavity is driven by the mechanical 
energy MGw(i)V'Jw/2 deposited into the IGM by the wind, 
but does not depend separately on Mgw and Vow - Hence, to 
determine Mgw, we must make an additional assumption. 
There are two limiting cases: One limit consists of having 
Mgw constant, in which case increasing the number of SNe 
will increase the wind velocity. The opposite limit consists 



of having a constant wind velocity Vgw, in which case an 
increase in the number of SNe results in a larger amount of 
matter being ejected. 

It would take detailed high-resolution simulations to de- 
termine which of these limits is correct. Ultimately, the criti- 
cal factor should be the spatial distribution of SNe. A single 
SN will only affect the ISM located in its vicinity. In the 
case of several SNe, there collective effect should critically 
depend on their level of clustering. If all SNe are concen- 
trated in a same location, the same region will be affected, 
and the net effect will be to eject the same matter, but at a 
larger velocity. If instead the SNe are distributed through- 
out the galaxy, each SN will affect a different part of the 
ISM, and the net result will be to eject more material, but 
at the same velocity. This last case is the limit we adopt 
in this paper. Gonsequently, the reader should keep in mind 
that our estimates of the amount of material ejected is an 
upper limit. Under this assumption the rate of mass loss by 
the galactic wind is proportional to the luminosity. 



Mowit) oc L{t) , 



(22) 



To determine the constant of proportionality, we first 
integrate the functions MGw(i) and L(t) over the lifetime 
of the galaxy: 



Jo 

E'°' = / Ldt, 
Jo 



(23) 
(24) 



where Mq^ is the total mass ejected into the galactic wind, 
and is the total energy deposited in the ISM. We can 
rewrite equation H22p as 



Mawit) 



M'a 



-L{t). 



(25) 



The problem is that A/qw and are not known until 
the simulation is completed, and to perform the simulation, 
we need to know these quantities in advance in order to 
calculate Mgw at every timestep. To solve this problem, we 
replace Mq-^ and in equation (|25|) by approximations 
that can be calculated ab initio, before actually performing 
the simulation. 



3.3.1 Estimate of E^°\ 

The luminosity L{t) is calculated as the simulation proceeds, 
but we can estimate it as follows: first, as we shall see below, 
the contribution of stellar winds to the luminosity becomes 
negligible once the SNe phase starts. If we neglect stellar 
winds, and also neglect the weak dependence of the SN lu- 
minosity on the metallicity (see Fig. [T]), we can directly es- 
timate the luminosity from the star formation rate M, . We 
define an integrated mass loss rate F{t) using: 



F(t + to 



M»{t')dt' , 



(26) 



where fonsot is the time elapsed between the formation of 
the stellar population and the onset of the first SN, and 
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Figure 6. Solid curve: SNc luminosity versus time, for a 10^ Mq 
galaxy witii an exponential SFR and a star formation efficiency 
/, = 0.1. Dotted curve: Integrated mass loss rate F{t), for the 
same galaxy. 



Figure 7. Escape fraction versus redshift, for a IO^Mq galaxy 
with a star formation efficiency /* = 0.1. 



^active is the time duration of the SNe phase0 Figure [6] 
shows the luminosity L{t) obtain from the simulation, and 
the quantity F{t) calculated using ionsct ~ 3.11 Myr and 
tactivc = 39.2 Myr (Fig. [l}. To a very good approximation, 
F{t) is proportional to L{t). We can therefore approximate 
equation (|25|) as 



Mg\ 



where 



M, 



GW 



F{t), 



Fdt. 



(27) 



(28) 



Both F{t) and _F'°' are calculated at the beginning of the 
simulation. 



3.3.2 Estimate of M^^. 

We still need to determine the total mass ejected by the 
galactic wind, Mq^, to be able to use equation (|27[1 . Like 
the total energy deposited in the ISM, Mqw is not 

known until the simulation is completed. To estimate it, we 
replace Mow in equation (|23p by an approximation that can 
be calculated at the beginning of the simulation. We then 
integrate to get Mq^, we substitute that value in equa- 
tion (|27[) . which then provides the mass loss by galactic wind 
during the simulation. To find an initial approximation for 
Mow, we first notice that observations at different redshifts 
suggest a relation between t he mass loss by galactic winds 
and the star formation rate (|Martinlll999l '). often expressed 
in terms of the ratio 



Mow 



(29) 



This value appears to vary significantly among galax- 
ies, with values r anging from 0.01 to 10 l|Veilleux et al.l 
I2OO5I ). [Murray e"t~a l. (2005) derived analytical relations be- 
tween the factor 77 and the velocity dispersion ct, for both 
momentum-driven and energy-driven winds. We focus in this 



^ The lifetimes of the shortest-lived and longest-lived progenitors 
are therefore tonsct and tonsct -I- iactivc , respectively. 



paper on energy-driven winds, but will consider momentum- 
driven winds in fu ture work. For energy-driven winds, 
[Murray et al.l l|2005h derived the following relation: 



Mgw = M. Co.iEs 



300 kms" 



(30) 



where g is the velocity dispersion, £3 = 1000ii''°*/Af*c^, and 
^0.1 = /™/0.1, with fw{Mgai) the fraction of energy provided 
by stars t hat is used to power th e wind, for a galaxy of 
mass Mgai (jScannapieco et al]l2002[ ). To calculate £3, we use 
Starburst99 with a 10® M© stellar population and a standard 
IMF. The total energy produced by SNe and stellar 

winds is always of the order of 10^"" '^ ergs, for all metallicities. 
This gives £3 — 0.011. Equation pO[) reduces to 



Mgw = 0.11M,/^ 



300 kms~ 



(31) 



For a galaxy of mass A/gai, t he velocity dispersion a is ca l- 
culated using the equation of lOppenheimer fc Pavel (|2008h : 



200 



5 X 1012Mq 



Ho 



1/3 



kms 



(32) 



where Zgt is the formation redshift of the galaxy, H is the 
Hubble parameter, with Ho being its present value, and h = 
//o/lOOkms-^Mpc-^ Equations (EIJ and ^ give us an 
estimate of Mgw. We then apply equation (|23p to calculate 
^^GWi which we substitute in equation (|27|) . This equation is 
then used to calculate the mass loss by galactic wind during 
the simulation [eq. (|14p ]. 

The parameter /osc is defined as the fraction of the ISM 
mass that escapes the galaxy: 



/esc — 



^0 

Q,hO Mgal 



(33) 



Figure [7] shows /esc versus formation redshift z^f, for various 
galactic masses Mgai, with a star formation efficiency of /* = 
0.1. The dependence on z^i comes entirely from the factor 
H{zgf) in equation (|32|l . Galaxies that form earlier have a 
larger velocity dispersion a for a given mass Mgai. Lower- 
mass galaxies eject a larger fraction of their ISM than higher 
mass galaxies, and for a 10^ Mq galaxy, we get /esc > 1, 
which simply means that the entire ISM will be ejected from 
the galaxy (so the actual /esc is unity). 
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4 GALACTIC WIND MODEL 

4.1 The Dynamics of the Expansion 

Tegmark et al. (1993, hereafter TSE) presented a formula- 
tion of the expansion of isotropic galactic winds in an ex- 
panding universe. In this formulation, the injection of ther- 
mal energy produces an outflow of radius R, which consists 
of a dense shell of thickness Rd containing a cavity. A frac- 
tion 1 — /m of the mass of the gas is piled up in the shell, 
while a fraction fm of the gas is distributed inside the cavity. 
We normally assume 5 ^ 1, /,„ <C 1, that is, most of the 
gas is located inside a thin shell. This is called the thin-shell 
approximation. 

The evolution of the shell radius R expanding out of a 
halo of mass A/gai, is described by the following system of 
equations: 



R = 



8tvG{p - Pext) 

Q.bH'^R 



HRy 



nn^R 



7?2 



L 



2nR3 



5Rp 
R 



(34) 
(35) 



where a dot represents a time derivative, Q,, flb, and H are 
the total density parameter, baryon density parameter, and 
Hubble parameter at time t, respectively, L is the luminosity, 
p is the pressure inside the cavity resulting from this lumi- 
nosity, and pcxt is the external pressure of the IGM. The four 
terms in equation (|34p represent, from left to right, the driv- 
ing pressure of the outflow, the drag due to sweeping up the 
IGM and accelerating it from velocity HR to velocity 7?, and 
the gravitational deceleration caused by the expanding shell 
and by the halo itself. The two terms in equation (|35p repre- 
sent the increase in pressure caused by injection of thermal 
energy, and the drop in pressure caused by the expansion of 
the wind, respectively. 

The external pressure, Pcxt, depends upon the density 
and temperature of the IGM. As in PMG07, we will assume 
a photoheated IGM made of ionized hydrogen and singly- 
ionized helium (mean molecular mass ^ = 0.611), w ith a 
fixed temperature Tigm = 10*K (jMadau et al.ll200ll ) and 
an IGM density equal to the mean baryon density pb- The 
external pressure at redshift z is then given by: 



Pcxt(2) = 



(36) 



The luminosity L is the rate of energy deposition or 
dissipation within the wind and is given by: 

L(t) — /u,(I/SNo + Lsw) — -'^comp , (37) 

where Lsnc and Lsw are the total luminosity responsible for 
generating the wind, as given by equations (|13p and p2p . 
respectively. Lcomp represents the cooling due to Compton 
drag against CMB photons and is given by: 



Lcomp = — -r (1 + -Z) pR 

4:5 nie \ he J 



(38) 



where at is the Thomson cross section, and T~,o is the present 
CMB temperature. 

The expansion of the wind is initially driven by the 
luminosity. After the SNe turn off, the outflow enters the 



"post-SN phase. 'Ifl The pressure inside the wind keeps driv- 
ing the expansion, but this pressure drops since there is no 
energy input from SNe. Eventually, the pressure will drop 
down to the level of the external IGM pressure. At that 
point, the expansion of the wind will simply follow the Hub- 
ble flow. 



4.2 Metal Distribution inside the Galactic Wind 

In the TSE model, the baryon density inside the cavity is 
Pi = Pb{t) fm / — S)"^, while the baryon density inside the 
shell is ps = pb{t)(l — /m)/[l — (1 — 3)^]. This gives a mass 
M = A-kR? pb{t) /i inside the volume of radius R, which is 
precisely the mass of the IGM within that radius in the ab- 
sence of a wind. Therefore, in the TSE model, the material 
inside the shell is swept IGM material, while the material in- 
side the cavity is IGM material left behind. The mass Mqw 
added by the galactic wind is neglected in the TSE model. 
Hence, the TSE model does not predict the distribution of 
that mass inside the cavity. This means that any distribu- 
tion we chose would not violate the assumptions on which 
the TSE model is based. 

The simplest approximation for the distribution of met- 
als in the wind consists of assuming that the metals car- 
ried by the galactic wind are sp read eve nly inside the cav - 
ity (see lScannapieco et al.ll2002l : PMG07: lBarai et al.ll201ll ') . 
This poses a problem for the metals ejected near the end 
of the post-SN phase, just before the wind joins the Hubble 
flow. These metals would have to be carried across the en- 
tire radius of the cavity, at velocities that exceed the wind 
velocity. Processes such as turbulence and diffusion could 
homogenize the distribution of metals inside the cavity, but 
only over a flnite time period. In this paper, we take the 
opposite approach, by assuming no mixing. Hence, the gas 
that escapes the galaxy early on will travel larger distances 
than the gas that escapes later. Since the metallicity and 
composition of the ISM evolves with time, the galactic wind 
will acquire both a metallicity gradient and a composition 
gradient, with the inner parts containing a larger propor- 
tion of metals. To simulate such wind, we use a system of 
concentric spherical shells. At the end of every timestep, the 
code calculates the amount of gas that will be added to the 
galactic wind: 



AMGw(iO = MGw{ti)At, 



(39) 



where ti is the time corresponding to the timestep. After the 
flrst time step, the wind reaches a radius R\ = 7?Gw(At). 
We deposit the wind material produced during that time 
step into the sphere of radius i?Gw(At), which constitutes 
our central shell. After the second timestep, the wind now 
reaches radius R2 = Raw {2 At). We first transfer the wind 
material located between and _Ri into a shell of inner ra- 
dius _Ri and outer radius R2, and we then deposit the wind 
material produced during the second timestep into the cen- 
tral shell. This process is then repeated. At every timestep 
n, a new shell is created between radii R„ and Rn-i, all the 
wind material is shifted outward by one shell, and the new 
material is deposited into the central shell. Finally, when the 



The stellar winds are still on, but their contribution is negligible 
at this point. 
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Figure 8. Mass of some elements present in ttie ISM versus rcd- 
sliift, for a 10^ A/q galaxy with a star formation efficiency /, =0.1 
and an exponential SFR. The colors corresponds to various ele- 
ments, as indicated. Solid lines: simulation without galactic wind; 
dashed lines: simulation with galactic wind. 



Figure 9. Mass loss rate of stars versus redshift, for a IO^Mq 
galaxy with a star formation efficiency /* =0.1. The various col- 
ors represent different SFRs, as indicated. Solid lines: simulations 
with SNe only; dotted lines: simulations with SNe and stellar 
winds. 



wind enters the post-SN phase, we no longer deposit materiel 
into the wind, and the shells expand homologously with the 
cavity. One nice feature of this model is that it has no free 
parameter. In particular, if does not depend on the value of 
the timestep. Using a different timestep would change the 
resolution at which the wind profile is determined, but not 
the profile itself. 



5 RESULTS 

Here we use the algorithm described in §4 to study the evolu- 
tion of starburst galaxies, in a concordance ACDM universe 
with density parameter fio ~ 0.27, baryon density parame- 
ter Qbo = 0.044, cosmological constant Aq = 0.73, and Hub- 
ble constant Ho = 71 kms~^Mpc~^ (h = 0.71). Because the 
parameter space is large, we focus on a fiducial case: a dwarf 
galaxy of mass Mgai = 10^ Mq forming at 2gf = 15. This case 
is particularly important because the vast majority of galax- 
ies in the universe are dwarfs, and in CD M cosmology, these 
galax ies tend to form at high redshift (e.g. iBlumenthal et ahl 
1 19841 ). Also our model assumes that galaxies form by mono- 
lithic collapse and not by the merger of well-formed galax- 
ies, an assumption that is more appropriate for dwarfs (e.g. 
IBlumenthal et al]|l984h . Note that the value of Zgc matters 
in the model. It affects the expansion of the galactic wind, 
which in turns affects the evolution of the ISM. 

We performed two simulations, one with our basic 
model, and another one in which we turned off the galactic 
wind. Figure [5] shows the abundances of a few elements in 
the ISM. Most of the ISM enrichment occurs between red- 
shifts z = 15 and z = 13, during the epoch of intense SNe 
activity. At lower redshifts, the enrichment by stellar winds 
dominates. The effect of the galactic wind is very small. 
Adding the wind results in a 10% increase in ISM metal- 
licity, caused by the removal of low-metallicity ISM during 
the early stages of the wind. The effect the galactic wind 
can be much more significant but this requires a smaller 
galactic mass Mgai or a larger star formation efficiency /, , 
or SFR much more extended in time than the ones we have 
considered. 

In the next three subsections, we explore the parameter 
space by varying, respectively, the SFR, the star formation 
efficiency /*, and the mass Afgai of the galaxy. 



5.1 Star Formation Rate 

Figure[9]shows the mass returned to the ISM by stellar winds 
and SNe, versus redshift, for the different SFRs. With an in- 
stantaneous SFR, there is only one stellar population and 
the mass loss profile is identical to the one provided by SB99. 
Stellar winds are absent in this case because all the stars 
were formed in a metal-free ISM. For the constant SFR, the 
increase in mass returned to the ISM is caused by the for- 
mation of more and more stars. Since SNe dominate over 
stellar winds, a plateau is eventually reached when the time 
of the simulation is equal to the lifetime of the SNe for the 
first generation of stars. After that moment, the contribu- 
tion of a new population is compensated by the death of 
an old population. The processus is the same for the expo- 
nential SFR, except that the mass of the stellar populations 
decreases with time. Hence, the death of an old population is 
replaced by the birth of a less-massive population, which ex- 
plains the absence of a plateau. For the constant SFR, there 
is a sudden drop at 2 = 10.7 which corresponds to the last 
SNe explosions. The material ejected after that corresponds 
to the giant phase of low-mass stars. 

Figure [To] shows the luminosity, internal pressure, and 
comoving radius of the galactic wind, for the various SFRs. 
Again, stellar winds do not have much effect on the results. 
One interesting aspect is that an extended period of star for- 
mation tends to produce a larger final radius for the outflow, 
compared with an instantaneous SFR, even though the total 
stellar mass A/» formed is the same. The top panel of Fig- 
ure [TT] shows the density profiles of the galactic wind. The 
density gradients are very strong, with the density dropping 
by 3 — 5 orders of magnitude from the center to the edge. 
This is caused mostly by the dilution resulting from the ex- 
pansion. The outer density profile is lower for the constant 
SFR than for the instantaneous and exponential ones. In 
our galactic wind model, the outer parts of the wind con- 
tain gas that was expelled by the galaxy at early time. The 
amount of material ejected during the early phases will de- 
pend of the mass loss rate Mgw at that time, which is pro- 
portional to L{t) [eq. (|25p ]. As Figure [TOI shows, in the early 
phases, L{t) is larger for the instantaneous and exponential 
SFR's than for the constant SFR, which leads to a larger 
amount of material being ejected, material which ends up 
in the outer parts of the wind. The dashed line shows the 
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Figure 10. Luminosity (top panel), internal pressure (middle 
panel), and comoving radius of the galactic wind (bottom panel) 
versus redshift, for a 10^ Mq galaxy with a star formation effi- 
ciency /» =0.1. Colors and linetypes have the same meaning as 
in Figure[9] The solid black lines show, for comparison, the model 
of PMG07. The dotted black line in the middle panel shows the 
external pressure of the IGM. 



Figure 11. Density profile (top), cumulative mass profile (mid- 
dle), and metallicity profile (bottom) of material ejected into the 
IGM by the galactic wind, at z = 0, for a ICp Mq galaxy with a 
star formation efficiency /* =0.1. Colors and linetypes have the 
same meaning as in Figures [9] and 1101 The dashed curves show 
the results for an exponential SFR with minimum SNe progeni- 
tor masses Mi = GMq and Mi = IOMq. In the top panel, the 
horizontal dashed line shows the density of the IGM inside the 
cavity, assuming fm = 0.1. 



density of the IGM inside the cavity, assuming fm = 0.1. 
Not surprisingly, the wind density exceeds the IGM density 
inside the galaxy, or immediately outside it. But at larger 
radii, the wind density drops several orders of magnitude 
below the IGM density. We calculated the mass of the IGM 
inside the cavity, assuming a shell thickness 5 = 0.05. For 
the 5 cases plotted in Figure llll the values are in the range 
3.86 - 7.27 X 10* M0. The mass added by the galactic wind 
is 2.73 X IO^Mq, or between 3.8% and 7.1% of the mass in 
the cavity. This justifies a posteriori the assumption made 
by TSE that the mass added by the wind can be neglected. 

The middle panel of Figure [TT] shows the cumulative 
mass profile of the galactic wind, that is, the mass Mawir) 
between and r. Even though the density is maximum in the 
center, the actual amount of ejecta located near the galaxy 
is negligible. For a galaxy of mass IO^A^q, collapsing at red- 
shift z — 15, the virial radius is r2oo = 2kpc, and the radius 
of the stellar component is even smaller. Essentially all the 
gas contained in the galactic wind has been ejected from the 
galaxy, and most of it is located at radius r > 100 kpc. 

The bottom panel of Figure [11] shows the metallic- 
ity profiles. The metallicity gradients have a different ori- 
gin, since the dilution caused by the expansion of the wind 
equally affects metals, hydrogen, and helium. Since the ma- 
terial located in the outer parts of the wind was ejected ear- 
lier than material located in the inner parts, the metallicity 
gradient simply reflects the time-evolving chemical compo- 
sition of the ISM. With an instantaneous SFR, the ISM is 
enriched in metals very rapidly. Hence, the gas ejected into 
the galactic wind at early times is already metal-rich. As a 



result, the metallicity in the outer parts of the wind is larger 
for the instantaneous SFR than for the other SFRs. 

The dashed curves in Figure[TT]show the effect of chang- 
ing the minimum SNe progenitor mass (for an exponential 
SFR). Lowering the minimum mass from 8Mq to 6Mq in- 
crease the final radius of the outflow by 12% and the total 
mass ejected by 59%. Increasing the minimum mass to IOMq 
reduces the final radius of the outflow by 8% and the total 
mass ejected by 21%. 

5.2 Efficiency of the Star Formation 

Figure [12] shows the effect of varying the efficiency of star 
formation for a IO^Mq galaxy with an exponential SFR. 
Apart from the fact that the mass loss rate increases with 
the number of stars formed, the most striking feature of this 
figure is that, for /« ^ 0.5, the galactic wind can be suffi- 
ciently powerful to eject the totality of the ISM. Eventually, 
the ISM is replenished by SN ejecta and stellar winds pro- 
duced by stars already formed. Figure [131 shows that stellar 
winds are more significant with /, = 0.5 than with /, = 1. 
Lowering /« spreads star formation over a longer period of 
time, enabling stars to form in an environment richer in 
metals, and resulting in stronger stellar winds. 

Figure [14] shows the evolution of the ISM metallicity. 
The metallicity increases faster with a higher star formation 
efficiency, since there are more stars available to enrich the 
ISM. When all the gas in the galaxy is eventually ejected into 
the IGM, the metallicity experiences a sudden increase be- 
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Figure 12. Rate of mass loss by SNe (blue), stellar winds (red), and galactic wind (black) versus redshift, for a W^Mq galaxy with 
an exponential SFR. Solid lines: simulations with SNe only; dotted lines: simulations with SNe and stellar winds. The various panels 
correspond to different star formation efficiencies /», as indicated. In the top panels, the solid and dotted lines are indistinct for the 
galactic winds and SNe. 




Figure 13. Rate of mass converted to stars versus redshift, for a 
10^ Mq galaxy with an exponential SFR. Solid lines: simulations 
with SNe only; dotted lines: simulations with SNe and stellar 
winds. The various colors correspond to different star formation 
efficiencies /, , as indicated. For /« = 0.2 and 0.1, the solid and 
dotted lines are indistinct. 



Figure 14. Metallicity of the ISM versus redshift, for a 10^ Mq 
galaxy with an exponential SFR. Colors and linetypes have the 
same meaning as in Figure [TS] 



fore reaching a maximum value. This sudden increase occurs 
when the mass remaining into the ISM becomes similar to 
the mass returned by stars. Then, the mass of the ISM keeps 
dropping, and the metallicity approaches the value corre- 
sponding to the last stellar ejecta. Afterward, the metallic- 
ity decreases because the last SNe ejected fewer and fewer 
metals. Then, when stellar winds are taken into account, 
the metallicity of the ISM continue to decrease with time 
because low-mass stars in their giant phases eject material 
composed mostly of hydrogen and helium. 

Figure [15] shows the evolution of the composition of 
the ISM. The importance of stellar winds becomes naturally 
larger when the star formation efficiency increases. The stel- 
lar winds do not have much effect on the total mass of the 



ISM. However, there is a significant difference in the abun- 
dances of carbon and nitrogen for /, > 0.1. Figure [T6l shows 
the density of various elements inside the galactic wind. The 
external part of the galactic wind has the composition of the 
ISM during the early phases. To have a significant effect on 
the composition of the external regions of the galactic wind, 
which are the prime contributor to the IGM enrichment, 
the enrichment of the ISM must happen rapidly, which is 
the case when /, is large. Figure [17] shows th metallicity 
profile of the galactic wind, for the various values of /» . The 
amount of metal ejected into the IGM seems large enough 
to fit observations. Various studies indicate that the IGM 
metallicity at redshifts between 2.5 and 3. 5 is of the order 
of 10~^ '^^a (iSorigaila fc Cowid [T996. : .HeUsten et all ll998l : 
iRauch et allll997l : lDave et al.lll99^ 
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Figure 15. Mass of various elements present in tiie ISM versus redshift, for a 10'*Mq galaxy witti an exponential SFR. The various 
panels correspond to different star formation efficiencies /, , as indicated. The colors corresponds to various elements, as indicated. Solid 
lines: simulations with SNe only; dotted lines: simulations with SNe and stellar winds. 
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Figure 16. Density profile of material ejected in the IGM, for a IQPMq galaxy with an exponential SFR. Black lines: total density; 
colored lines: density of various elements, as indicated. Linetypes have the same meaning as in Figure [151 
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5.3 The Mass of the Galaxy 

Usually, the evolution of the ISM is unafTected by the mass of 
the host galaxy. Increasing the mass of the galaxy increases 
the mass of the ISM, the mass in stars, the amount of gas 
ejected by SNe and stellar winds, and the amount of met- 
als ejected by exactly the same factor. The only thing that 
might affect this tendency is the mass loss caused by the 
galactic wind. If the power of the wind is moderate, the evo- 
lution of the ISM will be unaffected. This is the case for the 
most massive galaxies, because the energy deposited is less 
and less coherent, which reduces the fraction of energy /„, 
used to produce the galactic wind. For less massive galaxies, 
the ISM is more enriched, because the galactic wind expels 
more gas, which increases the relative importance of metals 
returned by stars. 

Figure [T5] shows the comoving radius of the galactic 
wind, for galaxies of various masses, all having /, = 0.1 
and an exponential SFR. The final radius increases with 
the mass, but this effect is weak and gets weaker at larger 
masses. The radius of the galactic wind R increases by a fac- 
tor of 1.6 from 10*^1/0 to 10® Mq, and 1.13 from 10® M© to 
10^° Mq. This results from the competition between several 
effects. The energy deposited into the ISM by SNe and stel- 
lar winds increases linearly with i\fgai, but the fraction /,„ of 
that energy which is used to power the wind decreases with 
Mgai. At large masses, oc 1/Mgai and the two effects can- 
cel out. At smaller masses, /m decreases slower than 1/Mgai 
and the energy available to power the wind increases with 
mass. That energy must compete with the gravitational pull 
of the galaxy [last term in eq. (|34|l ]. which increases with 
galactic mass. This effect reduces further the final radius 
of the wind, and at large masses R actually decreases with 
increasing mass. We do not include masses Mgai — 10^^ M© 
and 10^^ Mq in Figure fTSl because the galactic wind does not 
even start for those objects. With an instantaneous SFR, 
we can maximize the effect of the energy deposition and 
produce a wind from Mgai = IO^'^Mq, but this wind re- 
mains gravitationally bound to the galaxy and eventually 
falls back. 



6 SUMMARY AND CONCLUSION 

We have combined a population synthesis code, interpola- 
tion tables for the mass and composition of SN ejecta, and 
an analytical model for galactic winds into a single algo- 
rithm that self-consistently describes the evolution of star- 
burst galaxies. This model describes the evolution of the 
stellar populations in the galaxy, the evolution of the mass 
and chemical composition of the ISM, the propagation of 
the galactic wind, and the distribution and abundances of 
metals inside the galactic wind. In particular, the algorithm 
(1) provides a detailed calculation of the energy deposited 
into the ISM by SNe and stellar winds, which is respon- 
sible for driving the galactic wind, (2) takes into account 
the time-evolution of the chemical composition of the ISM, 
which directly affect the composition of the galactic wind, 
and (3) takes into account the removal of the ISM by galac- 
tic winds, which affects the metallicity of the ISM, and the 
metallicity of the stellar populations to follow. 

Our first results concern the SFR for the galaxy. For 
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Figure 17. Metallicity profile of material ejected in the IGM at 
z = 0, for a IO^Mq galaxy with an exponential SFR. The various 
colors represent difEerent star formation efficiencies . Colors and 
linotypes have the same meaning as in Figures 1131 and 1141 
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galaxies with an exponential SFR and a star formation efficiency 
/t = 0.1, for various galaxy masses M^^i in solar masses. Solid 
lines: simulations with SNe only; dotted lines: simulations with 
SNe and stellar winds. 



a given star formation efficiency /», a longer SFR tends to 
produce a galactic wind that reaches a larger extent, but 
this wind will be less dense. By increasing the star forma- 
tion efficiency, we can produce a wind that reaches a larger 
extent and has a higher metallicity near its front. In some 
cases, the energy deposited by the stars is sufficient to com- 
pletely expel the ISM. When it happens, star formation is 
shut down, and the galactic wind enters the post-SN phase 
prematurely. Hence, paradoxically, an increase in the star 
formation rate can sometimes result in a galactic wind that 
reaches a smaller extend. This happens with galaxies of 
masses Mgai = 10* or less, because their shallow po- 
tential well enables the complete removal of the ISM by the 
galactic wind. 

For galaxies with mass above 10^^ M©, the material 
ejected in the IGM always falls back onto the galaxy, no 
matter the value of /, . Therefore, in the case of energy- 
driven galactic winds, lower-mass galaxies are more likely 
to be the ones responsible for enriching the IGM and po- 
tentially perturbing the formation of nearby galaxies. Below 
10^^ M©, the extent of the galactic wind and its mass and 
metal content both increase with the mass of the galaxy 
at constant /*. With different values of /*, a less massive 
galaxy can sometimes produce a larger wind. 

Our current model does not take into account the ef- 
fect of Type la SNe. These are difficult to include, because 
of the uncertainties on the lifetime of the progenitors. The 
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simulations presented in this paper start at redshift z = 15, 
and end between redshifts 2 = 9 and 6. The corresponding 
time periods are shorter than 1 Gyr, which is shorter than 
the lifetime of several Type la progenitors. The energy pro- 
duced by Type la SNe is about 2 0% of the ener gy produced 
by Type II SNe (see Fig. 10 of iBensonI [20101 ). Hence, in- 
cluding the Type la SNe would result in a slightly larger 
final radius for the outflow. A Type la SNe can produce 
up to 7 times more iron t han a Type II SNe (see model 
W7 in iNomoto et"aLlll997l ). and their contribution to the 
iron enrichment of the ISM become important after 1 Gyr 
l|Wiersmall2O10l ). Hence, the abundances of iron we present 
in this paper are underestimated. But because of the delay, 
the additional iron produced would remain in the inner parts 
of the galactic wind. 

We have assumed a minimum value of Mi — 8Mq for 
the minimum mass of SNe progenitors. However, the correct 
value is actually quite uncertain. We did a few simulations 
with minimum masses of 6Mq and IOMq . Our preliminary 
results show differences of order 10% in the final radius of 
the outflow, and of order 20-60% in the total mass ejected, 
with the largest effect occurring when Mi is reduced. We 
intend to study this in more detail in the future. 

To conclude, properties of galactic winds depend on the 
host galaxy properties, such as the mass or star formation 
efflciency. The history of the ISM enrichment plays a de- 
terminant role in the chemical composition and extent of 
the galactic wind, and therefore its ability to enrich the 
IGM. The next step will consist of implementing this galac- 
tic outflow model into large-scale cosmological simulations 
of galaxy formation and the evolution of the IGM. These 
will be the first simulation of this kind to include a detailed 
treatment of the stellar winds and their impact on the chem- 
ical enrichment of the IGM 
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